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We have used density-functional-theory methods and the ab initio random structure searching 
(AIRSS) approach to predict stable structures and stoichiometries of mixtures of iron and oxygen at 
high pressures. Searching was performed for 12 different stoichiometries at pressures of 100, 350 and 
500 GPa, which involved relaxing more than 32,000 structures. We find that Fe 203 and FeC >2 are the 
only phases stable to decomposition at 100 GPa, while at 350 and 500 GPa several stoichiometries 
are found to be stable or very nearly stable. We report a new structure of Fe 2 03 with P2\2\2\ 
symmetry which is found to be more stable than the known RI 1203 (H) phase at pressures above 
~233 GPa. We also report two new structures of FeO, with Pnma and R3m symmetries, which 
are found to be stable within the ranges 195-285 GPa and 285-500 GPa, respectively, and two new 
structures of Fes 04 with Pca2\ and P2\/c symmetries, which are found to be stable within the 
ranges 100-340 GPa and 340-500 GPa, respectively. Finally, we report two new structures of Fe 40 s 
with P^/n and PSml symmetries, which are found to be stable within the ranges 100-231 GPa 
and 231-500 GPa, respectively. Our new structures of Fes 04 and Fe 40 s are found to have lower 
enthalpies than their known structures within their respective stable pressure ranges. 

PACS numbers: 91.60.Gf,64.30.-t,61.50.Lt,64.70.K- 


I. INTRODUCTION 

Fe/O is an important planetary mineral and a signifi¬ 
cant component of the Earth’s mantle, and it may also be 
important within the Earth’s core. Seismological studies 
have shown that the density of the Earth’s solid inner 
core is 2-3% too low to be plausibly explained by pure 
Fe, leading to the suggestion that light elements must 
also be present. 1 n IM! The identities and concentrations 
of light elements present in the core have yet to be es¬ 
tablished, although various arguments favour O, Si, C, S 
and H as the most likely candidatesP 

In an earlier paper we investigated the formation of 
iron carbides at high pressuresP and here we study 
iron oxides. O has been suggested as an important 
light element for explaining both the density jump at 
the Earth’s inner core boundary and the observed seis¬ 
mic wave velocitiesP The possible presence of O within 
the core has stimulated a range of experimental® and 
theoretical 9 investigations of the Fe/O system, albeit 
only for a handful of stoichiometries. In this paper 
we present a broader investigation spanning 12 different 
stoichiometries at three different pressures. We search 
for stable structures at each stoichiometry and pressure, 
and construct convex hull diagrams to determine which 
phases are stable against decomposition. We also present 
results for transitions between solid phases up to pres¬ 
sures of 500 GPa. 

Temperatures within the Earth extend to above 6000 
kP. 

Atoms such as O could be present in the Earth’s 
core as isolated impurities, as in a solid solution or, al¬ 


ternatively, Fe/O compounds could be formed. A solid 
solution has a much higher configurational entropy than 
an ordered Fe/O compound plus crystalline Fe, but the 
formation of crystalline Fe/O compounds could lead to 
more favorable bonding arrangements than are obtained 
in a solid solution. To help understand the conditions un¬ 
der which solid solutions and ordered compounds might 
be stable it is necessary to determine which compounds 
are energetically favorable, taking into account a wide 
range of stoichiometries. This is a formidable task in its 
own right, although it should be augmented by calcu¬ 
lating finite temperature effects arising from the vibra¬ 
tional motion of the atoms, the configurational entropy 
of the solid solution, and the electronic entropy. We leave 
the estimation of finite temperature effects to subsequent 
work. 

We have performed searches at 100, 350 and 500 GPa. 
The pressure at the bottom of the Earth’s mantle is about 
136 GPa, and our data at 100 GPa correspond to a pres¬ 
sure within the lower mantle. The pressures within the 
solid inner core are believed to be in the range 330-360 
GPajnmn and our results at 350 GPa are therefore rep¬ 
resentative of inner core pressures. The highest pressure 
at which we have performed searches of 500 GPa is above 
that found within the Earth, although such pressures oc¬ 
cur in other planets, including exoplanets. 

II. COMPUTATIONAL DETAILS 

Determining the most stable structure of a material 
corresponds to finding the global minimum energy state 
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within a high-dimensional space. A number of meth¬ 
ods have been developed to tackle this problem, which 
have met with varying degrees of success. The method 
we use, which has been shown to be successful in deter¬ 
mining structures that have subsequently been verified 
by experiment is called ab initio random structure 

searching (AIRSS)! 19 AIRSS has been applied to many 
systems at high pressures, including iron 20 'and oxygen. 21 
Within AIRSS the energy landscape is sampled by gener¬ 
ating random structures. The structures are then relaxed 
to the local enthalpy minimum for some chosen pressure 
using the castep plane-wave density functional theory 
(DFT) code. 22 This process is repeated until the lowest 
enthalpy structure is found several times or the available 
computational resources have been exhausted. 

To improve the efficiency of the searches we can impose 
a number of biases. One bias we use is to reject initial 
structures in which the smallest inter-atomic separation 
is less than some chosen value. Structures which contain 
atoms that are very close together are extremely high 
in enthalpy and may even cause the geometry optimiza¬ 
tion procedure to fail. Since the number of local min¬ 
ima is known to increase exponentially with the system 
size^, an additional bias may be required for searches 
over larger unit cells. We use a bias that exploits the 
fact that low enthalpy structures tend to possess symme¬ 
try. We can create an initial structure with a particular 
symmetry by first selecting a “structural unit” containing 
A formula units (fu). This structural unit is chosen to be 
the lowest enthalpy structure found in an A fu search. A 
random space group containing B symmetry operations 
is then chosen and applied to the structural unit. This 
generates a larger initial structure containing A x B fu. 
In this investigation A and B range from 1-3 and 2-4, 
respectively. 

Structures found to be metastable at one pressure 
may become stable at another. Expanding the enthalpy 
around the pressure po at which a search is performed 
gives 


dH 

H{p) = H(p 0 ) + (p-po) — 


P 0 


1/ , 2 d 2 H 

■ 2 (p - f, “ ) W 


( 1 ) 


P 0 


The second-order derivative in Eq. 0 is related to the 
bulk modulus and is computationally expensive to eval¬ 
uate. A useful approximation can, however, often be ob¬ 
tained by neglecting the second-order term and using the 
fact that the first derivative of H with respect to pressure 
p at po is equal to the volume vq, so that 


H(p) « H(p 0 ) +v 0 (p- po). (2) 

Since H(po) and Vo are evaluated within the structure 
searching, Eq. § may be used to estimate the enthalpies 
of the phases at any pressure p. This allows us to deter¬ 
mine approximate pressures at which phase transitions 
may occur, which we then refine by performing higher 
accuracy DFT calculations. 

We used ultrasoft pseudopotentiald^ to represent the 
cores of the Fe and O atoms. For the structure searching 


we used a fairly soft Fe pseudopotential in which only the 
4 s and 3d electrons were treated explicitly, while for O we 
treated the 2s and 2 p electrons explicitly. For the higher 
quality calculations we used a harder Fe pseudopotential 
in which the 3s, 3p, 4s, and 3d electrons were treated 
explicitly. Both the harder Fe pseudopotential and the 
O pseudopotential have been successfully tested up to 
terapascal pressures! 20 * 21 * 25 * 

The final structures were obtained using a two-step 
procedure. We first used AIRSS in medium-quality cal¬ 
culations which were optimized for computational speed. 
The medium quality settings consisted of a /c-point sam¬ 
pling grid of spacing of 27rx0.07 A -1 and a plane-wave 
cut-off energy of 490 eV. The lowest enthalpy structures 
for each stoichiometry were then further relaxed in a 
higher-quality calculation using a finer fc-point sampling 
grid of spacing 27rx0.03 A -1 and a plane-wave cut-off 
energy of 1000 eV. In all of our searches we used the 
Perdew-Burke-Ernzerhof (PBE) generalized gradient ap- 
proximation (GGA) exchange correlation functional.*- 2 ^ 

The Fe/O system is known to be problematic for DFT- 
based methods because of magnetic and strong correla¬ 
tion effects. The presence of magnetism in FeO at high 
pressures has been documented, although there are some 
conflicting reports. Mossbauer spectroscopy has shown 
that a continuous transition from high-spin to low-spin 
stated occurs at 90-120 GPa, leading to an eventual 
collapse at 140 GPa and 300 K. This is in contrast to ex¬ 
periments performed using X-ray emission spectroscopy 
for FeO at 300 K, which have suggested that a mag¬ 
netic state exists up to at least 143 GPsP^, while a mag¬ 
netic collapse in Fe 203 has been observed at 50 GPa™ 
In addition, zero-temperature static lattice GGA DFT 
calculations have predicted a magnetic collapse of FeO 
at around 200 GPa.® 

The existence of magnetic order under core conditions 
can be safely ruled out, as it would be destroyed by the 
high temperatures. We have not performed searches with 
spin-polarization because they are very costly, but to 
investigate the effects of magnetism on the enthalpies, 
we have re-relaxed our structures with high quality set¬ 
tings starting from a high spin state. If a relaxation did 
not achieve convergence we reduced the initial spin state 
and re-relaxed until convergence was achieved. We found 
that the small magnetic moments have a negligible effect 
on the enthalpies at the high pressures considered here 
(^0.002 eV per atom at 100 GPa in the worst case). For 
these reasons we do not believe magnetism plays a sig¬ 
nificant role at the pressures we have investigated and 
hence we have not considered it further. 

Similar arguments can be made about the effects of 
strong Coulomb correlations which arise in iron oxides 
from interactions involving the localized Fe d electrons. 
The strength of such correlations is conventionally mea¬ 
sured by the parameter U/W, where U is the Hubbard 
U energy and W is the d band width. Electrons tend 
to delocalize at high pressures, which has the effect of 
increasing W and reducing 17™ Under these conditions, 
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FIG. 1: Convex hull diagrams for Fe/O at 100, 350 and 
500 GPa. 


where the correlations are less strong, GGA DFT is likely 
to be reasonably accurate.^ 


III. CONVEX HULLS 


We have plotted our energy data on convex hull dia¬ 
grams to assess the stability of the phases with respect to 
decomposition. To construct the hulls we also need the 
enthalpies of the end members, pure Fe and pure O. Fe 
is predicted to adopt the hexagonal close packed (hep) 



X 

N 

*Sioo 

* 5*350 

* 5*500 

FegO 

0.10 

3 

[Amm2 (1)] 

[Amm2 (1)] 

[Amm2 (1)] 

FesO 

0.25 

4 

[P6m2 (1)] 

[P6m2 (1)] 

[PA/nmm (2)] 

Fe 2 0 

0.33 

6 

P&z/mmc (2) 

74/mmm (1) 

IA/mmm (1) 

Fe30 2 

0.40 

4 

Pm (3) 

C2/m (1) 

[Pm (2)] 

FeO 

0.50 

8 

P3 2 21 (3) 

R3m (1) 

R3m (1) 

Fe 4 0 5 

0.56 

2 

P4 2 /n (2) 

P3ml (1) 

P3ml (1) 

Fe30 4 

0.57 

4 

Pca2i (4) 

P2i/c (4) 

P2i/c (4) 

Fe 2 03 

0.60 

4 

Pbcn (4) 

P2i2i2i (4) 

P2i2i2i (4) 

Fe0 2 

0.66 

5 

Pa3 (4) 

Pa3 (4) 

R3m (1) 

FesOr 

0.70 

3 

P2/m (1) 

743d (2) 

743d (2) 

Fe03 

0.75 

4 

PI (4) 

PI (4) 

Cmcm (2) 

Fe0 4 

0.80 

4 

Fdd2 (2) 

P2i/c (2) 

P2i/c (2) 


TABLE I: The lowest enthalpy structures found from 
searching at 100, 350 and 500 GPa. X is the fraction of 
oxygen atoms. Entries which show space groups within 
square brackets indicate mixtures of structures of phases 
with other stoichiometries, see text. N is the maximum 
number of formula units used in the searches. Sp is the 
space group of the lowest enthalpy structure found at 
pressure P. In each case, the number of formula units 
in the primitive cell is given in round brackets. 


structure over the pressure range of interest here, while 
for oxygen at 500 GPa we use the P6 s/mme phase pre¬ 
dicted by Sun et al . 21 At pressures below 500 GPa we 
use the ( phased 

The enthalpy of formation per atom of a compound 
with respect to its elements is 


A H s = 


H s -(N Fe H Fe + N 0 H 0 ) 
No + TVpe 


( 3 ) 


where N x and H x are the number of atoms and enthalpy 
per atom of element x, respectively, and H s is the en¬ 
thalpy of the structure. A negative A H s indicates that 
the structure is stable with respect to decomposition into 
its elements, although it may be unstable with respect to 
decomposition into compounds of other stoichiometries. 
Structures lying above the convex hull are unstable with 
respect to decomposition into other compounds and are 
therefore less likely to form. Structures lying on the con¬ 
vex hull are thermodynamically stable and are therefore 
more likely to form. 

The convex hull diagrams at 100, 350 and 500 GPa 
shown in Fig. [l] exhibit fairly smooth variations with com¬ 
position. Figure Q|a)] shows that the only stoichiometries 
that lie on the hull at 100 GPa are Fe 2 03 and Fe0 2 . 
At 350 GPa, FegO, FesO, Fe 2 0, Fes0 2 , FeO, Fe 2 03 and 
Fe0 2 are found to lie on, or very close to the hull, and 
similar results are found at 500 GPa. Perhaps the most 
notable difference at 500 GPa is that Fe 04 is found to 
lie on the hull. Fe0 2 and Fe 2 03 are predicted to be sta¬ 
ble with respect to decomposition at all pressures investi¬ 
gated. Fe 304 , which is believed to be an important phase 
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(a) Fe20 (/4/mmra) (350 GPa) 



(b) FeO (RSm) (350 GPa) 



(c) Fe 2 0 3 (P2\2\2i) (350 GPa) (d) Fe0 2 (PaZ) (100 GPa) 



(e) Fe0 2 (mm) (500 GPa) 



(f) Fe0 4 (P2 1 /c) (500 GPa) 


FIG. 2: (Color online). Structures found from searching that lie on the convex hull. The O atoms are shown in red 
and the Fe atoms are in blue. 
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in the Earth’s lower mantle^, is found to be marginally 
unstable with respect to decomposition at the pressures 
considered here. 

Increasing pressure appears to stabilize Fe rich phases, 
an effect that is most noticeable between 100 and 350 
GPa. There is also a stabilization of O rich phases with 
increasing pressure. For example, FeC >3 and FeC >4 can 
be seen to move considerably closer to the hull between 
100 GPa and 500 GPa. The former also coincides with a 
transition from a low symmetry PI phase to a phase of 
Cmcmn symmetry. 

It is possible that structures that lie just above the 
convex hull correspond to metastable mixtures of phases 
rather than single phases. This is particularly likely for 
large unit cells. We have found several such cases in our 
results, and the fact that they lie just above the hull 
suggests that the interfacial energies between the differ¬ 
ent structures present are small. The P6m2 structure of 
FesO is found to be a mixture of phases consisting of Fe 
and Fe 20 while the PA/nmm structure is found to con¬ 
sist of Fe, Fe 2 0 and FeO phases. Finally, Fe g O (Amm2) 
is found to consist of Fe and Fe 2 0 phases and the Pm 
phase of Fes0 2 (at 500 GPa) is found to consist of Fe 2 0 
and FeO phases. These mixtures of phases are indicated 
by square brackets in Table [TJ 


IV. PHASE TRANSITIONS 


We predict pressure-driven phase transitions for all of 
the stoichiometries considered in this study, with the ex¬ 
ception of FegO. We have obtained accurate transition 
pressures by performing additional calculations at pres¬ 
sures around the values predicted using Eq. 0 - Plots 
of relative the enthalpies as a function of pressure are 
shown in Fig. [4] Three of the stoichiometries considered, 
Fe 2 03 , FeO and Fe 304 , have already been investigated 
extensively and their properties, particularly at ambient 
pressure, are well known. However, most of the stoi¬ 
chiometries that we discuss here have not been studied 
before and hence little is known about their high pres¬ 
sure behaviour. Unless stated otherwise, all of the lowest 
enthalpy phases were found by searching and, to the best 
of our knowledge, have not been discovered before. 

FegO. The most stable structure of FegO found has 
Amm2 symmetry, but it turns out to be a mixture of 
phases. 

Fe^O. See Fig. ^[a) Between 100 and 350 GPa, the 
most stable structure of FesO of P6m2 symmetry turns 
out to correspond to a mixture of phases rather than 
a single phase. At about 350 GPa we find a phase of 
P4/nmm symmetry to be the most stable, but it also 
turns out to be a mixture of phases. We also relaxed 
FesO between 100 and 500 GPa in the BiI 3 structure, 
which has previously been suggested as a low enthalpy 
phaseP However, we found this structure to be between 
0.3 eV and 2 eV per fu higher in enthalpy than our best 
structures found from searching. 


Pe 2 0. See Fig. |^[b) We predict that Fe 2 0 forms a 
structure of P 63 / rarac symmetry at 100 GPa. At about 
180 GPa we find a transformation to a phase of P3ml 
symmetry, which is stable up to about 288 GPa, at which 
point we predict it to transform into a phase of H/mmm 
symmetry. 

Fe 3 0 2 . See Fig. ^[c)| At 100 GPa we predict that 
Fes0 2 will adopt a structure of Pm symmetry. At about 
127 GPa, we predict a phase transition to a structure 
of C2/m symmetry followed by a further transition to a 
phase of Pm symmetry at about 486 GPa. 

FeO. See Fig. Jd) Previous DFT calculations for FeO 
have shown the stability of a phase of P3 2 21 symmetry^ 
between 65 GPa and 190 GPa. We also find th is stru cture 
from searching and our results shown in Fig. |^|[d)| agree 
well with those previously reported. However, Oganov 
et al predict that this phase will transform to the dis¬ 
torted NaCl structure at 190 GPa. We instead find that 
the P3 2 21 structure transforms to a new phase of Pnma 
symmetry at about 195 GPa, which then undergoes a 
transformation to another new phase of R3m symmetry 
at about 285 GPa. In addition to this, Fig. 11 shows 
another new phase of Cmcm symmetry, which we pre¬ 
dict to be stable at pressures above 500 GPa. Our results 
suggest a somewhat different picture from that obtained 
from compression experiments up to 324 GPa 4 and at 
temperatures up to 4880 K. These differences may arise 
for a number of reasons. Firstly, experiments on FeO 
normally use a non-stoichiometric compound, Fei-^O, 
where x is typically about 0.05, whereas our calculations 
assume stoichiometric FeO. Secondly, we have performed 
static lattice calculations in contrast to real-world experi¬ 
ments which may yield phases stabilized by nuclear vibra¬ 
tional motion. Thirdly, kinetic barriers between struc¬ 
tures could also influence the structures that are observed 
in experiments. Finally, the use of an approximate den¬ 
sity functional will lead to further uncertainties in our 
results. It is worth remarking, however, that the PBE 
density functional, which obeys the uniform limit and 
gives a good account of the linear response of the elec¬ 
tron gas to an external potentialP® i s likely to become 
more accurate as the pressure is increased. 

A notable example of where our predictions for FeO 
differ from observations is at relatively low pressures 
where experiments (extrapolated to 0 K) have shown that 
FeO adopts the rhombohedrally distorted NaCl (rBl) 
structure, which is predicted to be stable up to about 
100 GPa.® Above 100 GPa a phase transition to the NiAs 
(B 8 ) structure is predicted. We find the inverse-B 8 (iB 8 ) 
structure to be more stable than B 8 between about 100 
GPa and 165 GPa. The stability of the iB 8 structure 
has been predicted theoretically before 31 although it 
has yet to be verified by experiment. Although our re¬ 
sults correctly predict the B 8 phase to be more stable 
than the rBl phase at high pressures, our best structures 
obtained from searching are lower in enthalpy than the 
previously known phases across the entire pressure range 
considered. We are not aware of any experimental evi- 
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FIG. 3: (Color online). Electronic densities of states for structures lying on the convex hull. Each plot shows the 
total density of states and its decomposition into contributions from Fe and O. The Fermi energies are at 0 eV. 
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FIG. 4: (Color online). Plots of relative enthalpy against pressure for various stoichiometries. The reference 
structures appear as a horizontal line at 0 eV. Dashed lines represent structures known prior to this work. 


dence suggesting that a P3221 phase of FeO is stable at 
100 GPa, but we believe that this discrepancy may be 
explained by our arguments in the previous paragraph, 
and we note that earlier work suggested that the pre¬ 
dicted stability of the P3221 phase may simply be an 


artefact of the PBE density functional.^ 

Fe^O^. See Fig. ^e) Compression experiments have 
identified a phase of Fe 40 s with Cmcm symmetry to 
be the most stable between 5 and 30 GPa. 33 Our results 
show this phase to be unstable between 100 and 500 GPa 






































































































with respect to the new structures found in our searches. 
At 100 GPa we predict Fe^s to adopt a structure with 
P^/n symmetry prior to a transition to a structure of 
P3ml symmetry at about 231 GPa, which we predict to 
be stable to at least 500 GPa. 

Fe%O 4 . See Fig. ^T) Under ambient conditions FesCU 
adopts an inverse spinel structure of Fd3m symmetry!® 
Between ^20 — 60 GPa, Fes 04 is observed to make a 
slow transformation to a new phase (somet imes denoted 
by h-Fes 04 ) with orthorhombic symmetry !®^- 7 The ex¬ 
act structure of h-Fes 04 has, however, been difficult to 
resolve!® In light of this, a number of candidate struc¬ 
tures have been proposed. Fei et al. proposed a structure 
of Pbcm symmetry® which, although it agreed well with 
the available X-ray diffraction data, was later found to 
be inconsistent with Mossbauer spectroscopy data!® It 
was later suggested that a structure of Cmcm symmetry 
(sometimes referred to using the non-standard Bbmm 
setting) is also consistent with the X-ray diffraction data 
and it currently seems to be favoured as the high pres¬ 
sure phase of Fes 04. 39 First principles calculations have 
shown the Cmcm phase to be more stable than the Pbcm 
phase at pressures > 65 GPaf 37 Our searches at 100 GPa 
found a new phase with Pca2i symmetry to be the lowest 
in enthalpy up to ~ 340 GPa, above which a phase tran¬ 
sition is found to occu r to another new phase of P2\jc 
symmetry. Figure^?)] shows that the Cmcm phase is en¬ 
ergetically unfavourable at the pressures considered. The 
Pca2i structure is found to be much more stable than the 
Cmcm phase, being ^ 1.4 eV per fu lower in enthalpy 
at 100 GPa. The fact that the Pca 2 \ phase is also of or¬ 
thorhombic symmetry does not preclude it from being a 
candidate for h-Fes 04 . Ascertaining whether the Pca 2 \ 
phase is indeed h-Fes 04 would require a more detailed 
investigation, which we leave to subsequent work. 

Pe 2 03 . See Fig. Jg) At ambient pressures the stable 
phase of Fe 203 is hematite, which is a corundum-type 
structure of R3c symmetry™ The high pressure phase of 
Fe 203 has been proposed to be either a Rh 203 (II)-type 
structure (of Pbcn symmetry) or a GaFeOs orthorhom¬ 
bic perovskite structure!® Room temperature compres¬ 
sion experiments have shown that hematite transforms to 
the Rh 2 03 (H)-type structure at about 50 GPaf 41 How¬ 
ever, compression experiments at higher temperatures 
(800-2500 K) predict a phase transition to occur between 
hematite and the high pressure phase at a much lower 
pressure of ~ 26 GPa™ Ono et al. have suggested that 
the observed transition at room te mpera ture is between 
hematite and metastable RI 1203 (H).® 42 The perovskite- 
type phase is suggested to be kinetically inhi bited and 
therefore only obtainable at high temperatures!®® This 
view is supported by hybrid-DFT calculations, which 
show that the RI 1203 (H) phase is metastable with respect 
to the perovskite phase at 50 GPa!® At about 60 GPa 
and temperatures > 1200 K the high pressure phase is ob¬ 
served to trans form into a Calr 03 -type (post-perovskite) 
structure ! 40 ! 43 ! Hybrid-DFT calculations have predicted 
the anti-ferromagnetically ordered Calr 03 type structure 


to be the most stable between about 46 and 90 GPa with 
the low-spin RI 1203 (H) phase gradually becoming stable 
as the pressure approaches 90 GPa!® It appears likely 
that this phase will become stable at around 100 GPa and 
henc e become consistent with our results shown in Fig. 
^[g) We found the (non-magnetic) RI 1203 (H) structure 


of space group Pbcn in our searches at 100 GPa. Our 
calculations find this structure to be more stable than 
both the Calr 03 and GaFeOs (non-magnetic) structures 
between 100 and about 233 GPa. We predict that the 
RI 1203 (H) phase transforms into a new phase of P2\2\2\ 
symmetry at 233 GPa. 

Fe 02 - See Fig. ^[h) At 100 GPa Fe 02 we predict a 
structure of Pa3 symmetry, see Fig. S^d), which is stable 
across a wide pressure range of about 100-465 GPa. In 
this structure the Fe atoms form a face centred cubic (fee) 
arrangement. At 465 GPa the Pa3 phase is predicted to 
transform into a phase of R3m symmetry. 

FesOj. See Fig. |^[i)| Between about 100 and 160 GPa 
we predict a structure of P2/m symmetry, and above 
160 GPa we expect it to transform into a phase of /43d 
symmetry. 

FeOs. See Fig. ^ j) We find a structure of PI symme¬ 
try to be the most stable above 100 GPa, and we predict 
that it will transform into a different stucture of PI sym¬ 
metry above 162 GPa. At about 395 GPa this structure 
is predicted to transform into a phase of Cmcm symme¬ 
try. 

FeO±. See Fig. 11 We find a structure of Fdd2 sym¬ 
metry to be the most stable at 100 GPa, which is pre¬ 
dicted to transform into a structure of P2\jc symmetry 
at about 108 GPa. 

A summary of the structures found is given in Table 
[I] and the details of the structures are reported in the 
Supplemental Material!® 


V. DENSITY OF ELECTRONIC STATES 


An analysis of the electronic densities of states of the 
phases was performed with the LINDOS code® which uses 
a linear extrapolation method . 46 47 A Gaussian smearing 
of width 0.3 eV was applied to the density of states in all 
cases. Figure [3] shows the density of states for all of the 
structures found to lie on the convex hull. In all cases, 
with the exception of Fe 04 at 500 GPa, a significant den¬ 
sity of states is observed at the Fermi energy, indicating 
that these structures are metallic. 

The Fe 3s and Fe 3 p states lie at, respectively, ^87 and 
^54 eV below the Fermi energy in all Fe-0 phases. We 
found the O 2s states to lie at ^20 eV below the Fermi 
energy and almost all of the occupied electronic density 
of states around the Fermi energy arises from the Fe d 
bands, except in the O-rich structures. 

































9 


VI. CONCLUSION 

We have used DFT methods and the AIRSS technique 
to identify high pressure structures in the Fe/O system 
at zero temperature across a wide range of stoichiome¬ 
tries. This is a crucial preparatory stage for modelling 
the Fe/O system at finite temperatures. By constructing 
the convex hull at different pressures we have determined 
energetically favourable stoichiometries, some of which 
have not been considered before. Our results broaden 
the range of possible stable stoichiometries considerably 
from those studied previously. 

We have found new structures of FeO and Fe 203 which 
are calculated to be stable against decomposition at 350 
and 500 GPa. Structures with stoichiometries Fe 203 and 
Fe 02 are found to be stable against decomposition at 
each of the three pressures studied. FegO, FesO, Fe20, 
Fe 302 , FeO, Fe 203 and Fe02 are found to lie on, or very 
close to the hull at both 350 GPa and 500 GPa, although 
some of these structures are mixtures of phases, see Table 
E At 500 GPa Fe 04 is also found to lie on the hull. 

We have found structures of Fe 304 and Fe 405 at high 
pressures that are more stable than those predicted pre¬ 
viously, although they are unstable to decomposition at 
the pressures studied here. 


Increasing pressure tends to stabilize Fe/O compounds 
with respect to the separated elements over the range 
100-500 GPa. Increasing pressure stabilizes Fe rich 
phases, an effect that is most noticeable between 100 and 
350 GPa. There is also a stabilization of O rich phases 
with increasing pressure. 

The energy reduction from compound formation in 
Fe/O at high pressures is much larger than in Fe/C. For 
example, at 350 GPa, the DFT results reported in Ref. 
[3] show that the minimum in the convex hull occurs at 
Fe 2 C with an enthalpy reduction, compared with sepa¬ 
rated Fe and C solids, of about 0.26 eV per atom, while 
in Fe/O we find an enthalpy reduction of about 1.76 eV 
per atom for Fe02- Our results show that Fe/O com¬ 
pounds are stable (or nearly stable) over a wider range 
of stoichiometries at 350 GPa than Fe/C compounds. 


VII. ACKNOWLEDGEMENTS 

This work was supported by the Engineering and Phys¬ 
ical Sciences Research Council U.K. (EPSRC-GB). Com¬ 
putational resources were provided by the Cambridge 
High Performance Computing Service. 


1 G. Steinle-Neumann, L. Stixrude, R. E. Cohen, and O. 
Giilseren, Nature 413 , 57-60 (2001). 

2 G. Fiquet, J. Badro, E. Gregoryanz, Y. Fei, and F. Occelli, 
Phys. Earth Planet. Int. 172 , 125 (2009). 

3 G. L. Weerasinghe, R. J. Needs, and C. J. Pickard, Phys. 
Rev. B 84 , 174110 (2011). 

4 H. Ozawa, F. Takahashi, K. Hirose, Y. Ohishi, and N. 
Hirao, Science 334, 792 (2011). 

5 H. Ozawa, K. Hirose, S. Tateno, N. Sata, and Y. Ohishi, 
Physics of the Earth and Planetary Interiors 179 , 157 
( 2010 ). 

6 R. A. Fischer, A. J. Campbell, G. A. Shofner, O. T. Lord, 
P. Dera, and V. B. Prakapenka, Earth and Planetary Sci¬ 
ence Letters 304 , 496 (2011). 

7 M. Murakami, K. Hirose, S. Ono, T. Tsuchiya, M. Isshiki, 
and T. Watanuki, Physics of the Earth and Planetary In¬ 
teriors 146 , 273 (2004). 

8 D. Alfe, G. D. Price, and M. J. Gillan, Phys. Earth Planet. 
Inter. 110 , 191 (1999). 

9 D. Alfe, G. D. Price, and M. J. Gillan, Geophys. Res. Lett. 
27 , 2417 (2000). 

10 S. Anzellini, A. Dewaele, M. Mezouar, P. Loubeyre, and 
G. Morard, Science 340 , 464 (2013). 

11 B. J. Wood, Earth Planet. Sci. Lett. 117 , 593 (1993). 

12 L. Stixrude, E. Wasserman, and R. E. Cohen, J. Geophys. 
Res. 102, 24729 (1997). 

13 D. Alfe, M. J. Gillan, and G. D. Price, Earth Planet. Sci. 
Lett. 195 , 91 (2002). 

14 M. J. Gillan, D. Alfe, J. Brodholt, L. Vocadlo, and G. D. 
Price, Rep. Prog. Phys. 69, 2365 (2006). 

15 C. J. Pickard and R. J. Needs, Phys. Rev. Lett. 97, 045504 
(2006). 


16 C. J. Pickard and R. J. Needs, Phys. Rev. B 76, 144114 
(2007). 

17 A. D. Fortes, E. Suard, M.-H. Lemee-Cailleau, C. J. 
Pickard, and R. J. Needs, J. Am. Chem. Soc. 131, 13508 
(2009). 

18 G. I. G. Griffiths, A. D. Fortes, C. J. Pickard, and R. J. 
Needs, J. Chem. Phys. 136 , 174512 (2012). 

19 C. J. Pickard and R. J. Needs, J. Phys.: Condens. Matter 
23 , 053201 (2011). 

20 C. J. Pickard and R. J. Needs, J. Phys.: Condensed Matter 
21, 452205 (2009). 

21 J. Sun, M. Martinez-Canales, D. D. Klug, C. J. Pickard, 
and R. J. Needs, Phys. Rev. Lett. 108, 045503 (2012). 

22 S. J. Clark, M. D. Segall, C. J. Pickard, P. J. Hasnip, M. 
I. J. Probert, K. Refson, and M. C. Payne, Z. Kristallogr. 
220 , 567 (2005). 

23 F. H. Stillinger, Phys. Rev. E 59 , 48 (1999). 

24 D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

25 J. Sun, D. D. Klug, C. J. Pickard, and R. J. Needs, Phys. 
Rev. Lett. 106 , 145502 (2011). 

26 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77 , 3865 (1996). 

27 M. P. Pasternak, R. D. Taylor, R. Jeanloz, X. Li, J. H. 
Nguyen, and C. A. McCammon, Phys. Rev. Lett. 79 , 5046 
(1997). 

28 J. Badro, V. V. Struzhkin, J. Shu, R. J. Hemley, and H.-k. 
Mao, Phys. Rev. Lett. 83, 4101 (1999). 

29 M. P. Pasternak, G. Kh. Rozenberg, G. Yu. Machavariani, 
O. Naaman, R. D. Taylor, and R. Jeanloz, Phys. Rev. Lett. 
82 , 4663 (1999). 

30 R. E. Cohen, I. I. Mazin, and D. G. Isaak, Science 275 , 
654 (1997). 



10 


31 Z. Fang, I. V. Solovyev, H. Sawada, and K. Terakura, 
Phys. Rev. B 59 , 762 (1999). 

32 Y. Ma, A. R. Oganov, and C. W. Glass, Phys. Rev. B 76, 
064101 (2007). 

33 B. Lavina, P. Dera, E. Kim, Y. Meng, R. T. Downs, P. F. 
Week, S. R. Sutton, and Y. Zhao, Proc. Natl. Acad. Sci. 
USA 108 , 17281 (2011). 

34 L. S. Dubrovinsky, N. A. Dubrovinskaia, C. McCammon, 
G. Kh. Rozenberg, R. Ahuja, J. M. Osorio-Guillen, V. 
Dmitriev, H.-P. Weber, T. Le Bihan, and B. Johansson, 
J. Phys.: Condens. Matter 15 , 7697 (2003). 

35 A. R. Oganov, Y. Ma, C. W. Glass, and M. Valle, Psi-k 
Newsletter 84, 142 (2007). 

36 Z. Fang, K. Terakura, H. Sawada, T. Miyazaki, and I. 
Solovyev, Phys. Rev. Lett. 81 , 1027 (1998). 

37 S. Ju, T.-Y. Cai, H.-S. Lu, and C.-D. Gong, J. Am. Chem. 
Soc. 134 , 13780-13786 (2012). 

38 Y. Fei, D. J. Frost, Ho-Kwang Mao, C. T. Prewitt, and D. 
Haiisermann, Am. Mineral. 84 203-206 (1999). 

39 Camilla Haavik, Svein Stolenl, Helmer Fjellvag, Michael 


Hanfland, and Daniel Hausermann, Am. Mineral. 85, 514- 
523 (2000). 

40 N. C. Wilson and S. P. Russo, Phys. Rev. B 79, 094113 
(2009). 

41 G. K. Rozenberg, L. S. Dubrovinsky, M. P. Pasternak, O. 
Naaman, T. L. Bihan, and R. Ahuia, Phys. Rev. B 65, 
064112 (2002). 

42 S. Ono, K. Funakoshi, Y. Ohishi, and E. Takahashi, J. 
Phys.: Condens. Matter 17, 269 (2005). 

43 S. Ono and Y. Ohishi, J. Phys. Chem. Solids 66, 1714 
(2005). 

44 See Supplemental Material at (link) for details of the struc¬ 
tures obtained in this study. 

45 A. J. Morris and C. J. Pickard, LINDOS version 1.3 User 
Manual, University College London (2010). 

46 C. J. Pickard and M. C. Payne, Phys. Rev. B 59, 4685 
(1999). 

47 C. J. Pickard and M. C. Payne, Phys. Rev. B 62, 4383 

( 2000 ). 



